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Abstract 

The experimental determination of the partial width r{H — > 77) of an 
intermediate mass Higgs is among the most important measurements at a future 
photon photon collider. Recently it was shown that large non-Sudakov as well as 
Sudakov double logarithmic (DL) corrections can be summed to all orders in the 
background process 77 ( = 0) — > qq. It was found that positivity and stability 
of the cross section was only restored at the four-loop level. One remaining large 
source of uncertainty stems from the fact that the scale of the strong coupling 
is unspecified within the double logarithmic approximation. In this paper we 
include the leading and next-to-leading order running coupling to all orders. We 
thus remove the inherent scale uncertainty of both the exact one-loop and all- 
orders DL result without encountering any Landau-pole singularities. The effect 
is significant and, for the non-Sudakov form factor, is found to correspond to an 
effective scale of roughly as{9m'^). 
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1 Introduction 



The model independent knowledge of the total Higgs width combined with respective 
branching ratios allows an experimental determination of various partial widths of the 
Higgs boson 0. It is therefore of fundamental importance for a detailed under- 
standing of the electroweak symmetry breaking mechanism. For an intermediate mass 
Higgs, i.e. with a mass below the W^W~ threshold (in the MSSM mass window), a 
photon linear collider (PLC) offers so far the only possibility of a direct and model 
independent measurement of Ttot- Using Compton backscattering ||^, ^, Q of initially 
polarized electrons and positrons the partial width T{H — > 77) can be determined 
with very good precision at a PLC 0. This quantity is of significant interest in its own 
right, since all charged massive particles contribute in the loop which might therefore 
be sensitive to new physics. With the respective branching ratio BR(if — > 77), de- 
termined in measurements at the LHC (and possibly also the NLC), the total width 
can be reconstructed. 

The dominant non-Higgs background for this energy regime is 77 (J^ = 0) — > qq 

2 

with q = b,c. While this background is suppressed by ^ at the Born level (unlike 
the Jz = ±2 channel), higher-order QCD radiative (Bremsstrahlung) corrections in 
principle remove this suppression 0. In addition, large virtual non-Sudakov double 
log(s/m^) logarithms (DL) are present which at one loop can even lead to a negative 
cross section P, 0]. In Ref. 0] we presented explicit three-loop results in the DL ap- 
proximation which revealed a factorization of hard (non-Sudakov) and soft (Sudakov) 
double logarithms for this process and led to the all-orders resummation in the form 
of a confluent hypergeometric function 

^vM+soft = ^Born (l + 2^2(1, 1; 2, | ^^) + 2 ^ ^^^(l, 1; 2, | ^^)' 




.^°s;^i2-^°s4ifj+^°s4if-' + y 



(1) 



where JF = —Cpf^ log^ ^ is the one-loop hard form factor. The Born cross section is 
given by 

dcos9^' /V . ; -^7 ^ (l-/52cos20)' ^ ^ 



where (3 = yl — ^^ denotes the quark velocity and Qq the charge of the quark with 
mass rrig. a = ^ is the fine structure constant and y/s the center-of-mass energy of 
the initial photons. 

In Ref. |]rT| it was pointed out that one needs to include at least four loops (at the 
cross section level) of the non-Sudakov hard logarithms in order to achieve positivity 
and stability. A major source of uncertainty remained in the scale choice of the QCD 
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coupling, two possible 'natural' choices — tts(s) and as(mg) — yielding very different 
numerical results. 

In this work we include the exact leading and next-to-leading order running cou- 
pling into the derivation of both the massive Sudakov and the novel hard (non-Sudakov) 
form factors, thus removing a major source of error in both the exact one-loop and 
all-orders DL results. Our technique is based on the observation that in a particular 
loop in a multi-loop ladder, the correct scale for the strong coupling is the character- 
istic transverse momentum flowing round the loop. Details are given in Section ^ In 
Section § we present numerical results including the full one-loop radiative corrections 
based on Refs. |^ and with the new renormalization group improved form factor. 
Section ^ contains our conclusions. 



2 Renormalization Group Improved Form Factors 

In this section we will include the QCD running coupling into the DL form factors given 
in Eq. A very important result incorporated in this expression for the virtual plus 
soft real Bremsstrahlung cross section is the factorization between Sudakov and non- 



Sudakov double logarithms [|10|. For our purposes here this means that it is sufficient 
to consider the corrections for the basic topologies shown in Fig. |l| separately. The 
complete RG improved result is then given by the same factorized structure.^] 

Within the DL approximation the scale at which to evaluate the QCD coupling is 
unspecified, and this is therefore a major source of uncertainty in the determination 
of the size of the bb, cc background for intermediate mass Higgs production. In the 
following we will derive the renormalization group effect of inserting a running coupling 
into each loop evaluation using the exact one- and two-loop solution of the /^-function. 
All higher-order RG-terms will then be suppressed by -. — is-. We begin with the more 

familiar case of topology A (see Fig. |l|). 



2.1 The Sudakov RG-Form Factor 



In the derivation of the leading logarithmic corrections in Ref. |^ we used the familiar 
Sudakov technique |T^, 0, |T3| of decomposing loop momenta into components along 

^In effect, introducing the running coupling softens the contributions ft'om the loops. It can- 
not promote sub-leading logarithms to the leading-logarithm level, and therefore cannot spoil the 
factorization. 
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Figure 1: The schematic one-loop soft (Sudakov, "A") and hard (non-Sudakov,"B,C") 
topologies contributing to the DL form factor. These graphs are obtained from the 
one-loop box diagram in the process 77(^2 = 0) — > qq. The blob denotes a hard 
momentum flowing through the omitted propagator relative to the soft momentum / 
or li in the DL phase space. For higher-order DL contributions, only corrections to 
these topologies need to be taken into account. 



external momenta, denoted by {a,P}, and those perpendicular to them, denoted by 
l±. For massless fermions, the effective scale for Sudakov double logarithms for the 
coupling at each loop is ^^(I^), as was shown in Refs. [jl6|, Q by direct comparison 
with explicit higher-order calculations. For massive fermions the effective scale is also 
given by 1^ = —l\ > as the dominant double logarithmic phase space is given by 

— <^ — ^ 1 [M (on a formal DL level, setting m = X yields the massless Sudakov 



form factor). We start by writing 

— 2 ,2 / ^ 2^\2 12 — 12 

l + /5o^logi^ + A(^) log^ 1 + clogi 

where Pq = — -^Tpnp, Pi = — ^CaTf^f — ^CfTfUf and for QCD we 

have Ca = 3, Cf = ^ and Tj? = | as usual. Up to two loops the massless /?-f unction is 
independent of the chosen renormalization scheme and is gauge invariant in minimally 



subtracted schemes to all orders [|12|. These features will also hold for the derived 
renormalization group improved form factors below. 

In case of the Sudakov DL's, the DL-phase space structure is very transparent in 
terms of the parameters. This is shown in Fig. ^ for the virtual and virtual 

plus soft real Bremsstrahlung phase space. We therefore use the on-shell condition 
^ix ~ -scuA, even though the running coupling will now depend on two integration 
variables. 

The DL result for the virtual one-loop Sudakov form factor can be written in terms 
of integrals over ai and Pi as: 

^It^ = /; ^ C —^{saiPi - X^)Q{ai - ^Pi)Q{Pi - ^ai) 

(4) 
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Figure 2: The schematic DL-phase space for the massive virtual (upper diagram) and 
the virtual plus soft real Sudakov form factor (lower diagram) . The sum is gluon mass 
independent for A < ^M^. . 




Figure 3: A schematic Feynman diagram leading to the Sudakov double logarithms 
in the process jjiJz = 0) — > qq with i gluon insertions. The blob denotes a hard 
momentum going through the omitted propagator in the DL-phase space. Crossed 
diagrams lead to a different ordering of the Sudakov variables with all resulting Ca 
terms canceling the DL-contributions from three gluon insertions |T^. The scale of the 
coupling = 1^ is indicated at the vertices and explicitly taken into account in this 
work. 



Here A is a fictitious gluon mass introduced to regulate the infra-red divergences in the 
real and virtual gluon integrals. For the inclusion of double logarithms from the real 
Bremsstrahlung contribution we introduce a cutoff for the energy integration according 
to X < L < L, i.e. in terms of Sudakov variables, ai + di < 
for the real DL Bremsstrahlung contribution: 

asCp 



Thus we find 



DL 



TT 
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(5) 



assuming only A < These results are of course well known ||T3| , [T^, ^ but now 

can be used to insert a running coupling into each loop integration. The DL result 
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is sufficient for this purpose as the RG-logarithms are sub-leading at the next order 
in Ug. At this point we only give the result for the sum J^^^ + 2Tg^ , as only this 
combination is relevant for a physical cross section. We emphasize at this point that in 
a adopting a certain jet definition, we need to make sure that it does not restrict the 
exponentiation of the energy cut dependent piece of the soft gluon matrix elemtent. 
In the Appendix we give the result for the renormalization group improved form 
factor and in addition the two-loop results, thus checking the requirement that all 
one-loop RG-form factors exponentiate. We find 

-pRG.r^-pRG _ J / r dPi_ dai 




l + clog((^-A)A^ 



+ log — log h log ^ log 



1 , as{m^)as{2lc^/s) . 

An expansion in ^^(m^) gives the double logarithmic result plus subleading terms in Po 
etc. The dependence on the cutoff Ic will vanish when the hard-gluon Bremsstrahlung 
contributions are added order by order. Fig. ^ compares the second term in the ex- 
ponential obtained from the renormalization group improved result in Eq. ^ with the 
explicit two-loop result given in Eq. ^ of the Appendix.^] We emphasize that the two 
loop result depends on two "running" scales , li^^ and l2j^. The figure clearly shows 
that our one-loop result exponentiates as expected. 



Since in this work we are only able to include the hard one-gluon matrix elements 
(the exact NNLO corrections are presently unknown), the higher-order terms are at 
this point of a more academic interest. Our two jet definition (see below) restricts 
the higher-order hard gluon radiation phase space sufficiently so that it is reasonable 
to neglect more than one gluon emission. As mentioned above. Fig. ^ shows, using 
explicit RG-improved two loop DL-results in the Appendix, that the result in Eq. || 
exponentiates as expected. We do have to consider the subleading terms already 

^The two-loop numerical result is obtained by Monte Carlo integration with 10^ evaluations for 
each of the fifty iterations. 
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Real+Soft Exp.: 
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Figure 4: The second term in the exponential obtained from the renormahzation group 
improved massive one-loop Sudakov form factor compared with the explicit two-loop 
results given in the appendix. The latter depends on a different scale in each loop inte- 
gration. The quark mass is kept fixed at rrify = 4.5 GeV and the values for the real gluon 
energy cutoff Ic is indicated in the figure. The result in Eq. ^ clearly exponentiates. 
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entering at the one-loop level as the hard Bremsstrahlung contribution will also contain 
subleading terms in l^. The complete result for the renormalization group improved 
massive Sudakov form factor is thus given by: 
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1 + c log ( 













— log log ^ log — . s 

as{s) c y/s as (^^) 



- , ^svi^)as{'2lcVs) 1 s s 
log — - — 7:;rzjv- + o log — + log - 1 + — , (7) 



assuming only ^ -C 1. Expanding in ^^(m^) gives the DL-Sudakov form factor in 
Eq. |I| together with subleading terms proportional to Pq and subsubleading terms 
proportional to We emphasize that the two- loop running coupling is included in 
Eq. to all orders and that all colhnear divergences are avoided due to the fact that 
we keep all non-homogeneous fermion mass terms. We next turn to the virtual hard 
DL corrections and investigate the RG effects for these contributions. 

2.2 The Non-Sudakov RG-Form Factor 

In the absence of exact two-loop QCD corrections to the process 77(^2 = 0) — > qq it 
might seem unclear how to perform a renormalization group improvement for the hard 
non-Sudakov DL's. We find, however, following similar arguments as in the case of the 
massive Sudakov DL's of the previous section, that the same effective scale 1^^ is also 
valid for the novel non-Sudakov logarithms which were summed to give the confluent 



hypergeometric function 2-F2 in Ref. [IQ]. One way to see this is to note that on a formal 
level, all insertions of gluons into the hard topology shown in Fig. ^ have the same 
structure as in the Sudakov case. Only the last fermion loop integration separates the 
two cases by effectively regularizing soft divergences with the quark mass. The strong 
coupling receives no renormalization from this last integration so that the scales of the 
couplings at each order are determined by the same renormalization group arguments 
as for the Sudakov case. In this context it is worth pointing out that our analysis here 
also thus removes the largest uncertainty in the exact one- loop calculation of Ref . , 
as the scale of in that work was unspecified at this order. 



From the exact next-to-leading order result for the running coupling in Eq. ^ it is 
clear that a formulation in terms of 1^ of the series leading to the non-Sudakov double 
logarithms is more adaptable to a renormalization group improvement. In Ref. |TD] 
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Figure 5: The schematic Feynman diagrams leading to the hard (non-Sudakov) double 
logarithms in the process 77(^2 = 0) — > qq with i + 1 gluon insertions. The blobs 
denote a hard momentum going through the omitted propagator in the DL-phase space. 
Crossed diagrams lead to a different ordering of the Sudakov variables and are correctly 

2 

accounted for by a factor of (i + 1)! at each order. The scale of the coupling = |j 
is indicated at the vertices and included explicitly in this work. The topology on the 
left-hand diagram is Abelian like, and the one on the right is non-Abelian beyond one 
loop. 



these contributions were derived by first integrating over all l\ with integrals over 
{a,j3} pairs remaining. For the topology depicted on the left in Fig. ^we find that it 
can be reformulated as follows: 



^.F.(u;2.|;i^) A ^ /' f nr(„)/'^?^ rf±i(^ 

2 2 ^^Jrr^ a (3 f-J-Q J a, ai+i Jp, A+i V 27r 



/ , 12- M — n rfn 



where J-" = — "^^^ log^ ^ denotes the hard one-loop form factor of Ref. In this 
result the product is set to one for n = 1, and contains nested integrals for n > 2 
with Iq = 1^ and ao = a. From this expression it is clear that an incorporation 
of the running coupling in Eq. ^ will not contain any Landau pole singularity as 
< if^ < l^. We now include the running of as according to Eq. |^ as follows. For 
each gluon insertion we have 

logajlog(l + c log^^^ 



If+lx 1 + c log ^ ^ ^ 



log a, log (9) 
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With 



we find 



(10) 



i-l I 



It is clear from tliis derivation tliat for i-gluon iterations we fiave 

(-1)^+^ ^a.K)Vlog^a^^g, Mlil _ ^^2) 



i\ \ c J a as(m2) 

Tfius we finally arrive at the complete renormalization group improved result for the 
hard non-Sudakov form factor corresponding to the left ( "Abelian" ) topology in Fig. ^ 



I dll fC^y^W -a,{m') y a.(Ii) Ij a.(li) 

£;^y™2l^l2^j c Jt + lV.^^ ^^^^ 



i=0 



The effect of the renormalization group improved hard form factor is shown in Fig. |^ 
for the case of the 6-quark. The effective scale of the coupling in the DL approximation]^ 
is roughly 9m^. For the series obtained in the non- Abelian topology on the right of 
Fig. ^ only the color factor changes, as was shown in Ref . . The correct substitution 
at the {i + l)-loop level is 

C'+'-^cJ^\ (14) 



2 

with an additional factor of 2 as this topology occurs twice in the process 77 — > qq. 
In summary, the complete virtual renormalization group improved hard form factor is 
given by 

^ijG ^ ^ r dll fCFV^W c^sim') y «.(li) Ij 1,^. , 

fro^-^ li V27rj 1^ c ) (z + 1)! ^ s a,(Ii) + 

^ fs dll C,C\ ( as{m') \ ' asjll) , ,+1 ll , , «.(m^) 
^ / 12 22m^m 1^ c J (^ + 1)! ^ s ^ a.(li) ' ^ ^ 



j=0 



and thus 



^RG 



{l+^rye.p{^i^ + 2^i^) (16) 



''By this we mean the scale which when used for in the non-RG-improved DL result reproduces 
numerically the RG-improved result derived here. 
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where the RG-improved massive Sudakov form factor is given in Eq. ^ Fig. ^ shows the 
size of the effect at the cross-section level for m = nib = 4.5 GeV. We take as{ml) = 0.2 
and use Eq. ^ with n^? = 5 to arrive at as{s). The effect is significant and lies almost in 
the middle of the 'theoretical' upper and lower limits for the non-Sudakov form factor 
(suppressing the Sudakov term for now) as indicated in the figure. 



In the following section we will now investigate the effect on the full cross section, 
including the Sudakov form factor as well as the exact one loop radiative corrections. 



3 Numerical Results 



In Ref. we presented numerical predictions for an (infra-red safe) two-jet bb cross 
section in e~^e~ collisions in the energy range y/s = 100 — 200 GeV. We used a modified 
Sterman- Weinberg cone definition. Thus, at leading order (i.e. 77 —>■ bb) all events 
obviously satisfy the two-6-jet requirement.^ This defines our 'leading order' (LO) 
cross section. At next-to-leading order (NLO) we can have virtual or real gluon 
emission. For the latter, an event is defined as two-6-jet like if the emitted gluon 

either I. has energy less than e^/s, with e <^ 1, 

or II. is within an angle 26 of the b or b, again with 5 <^ 1. 

we further subdivided region I according to whether the gluon energy is greater or less 
than the infrared cutoff Ic (< e). Adding the virtual gluon corrections to this latter 
(soft) contribution, to give asv, and calling the remaining hard gluon contribution o"h, 
we have 

o-2j = crsv(^c) + cru{lc, e, 5) . (17) 



In Ref. |TI| we evaluated each part of this cross section exactly to 0{as) and in addition 
we included the resummed hard non-Sudakov form factor in ctsv- This was necessary 
to yield a positive cross section. 

We can now update these predictions using the RG-improved expressions for the 
resummed form factors. Thus 

(Tsv (18) 



We apply an angular cut of |cos6', t| < 0.7 to ensure that both jets lie in the central region. 
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Figure 6: The effect of incorporating a running coupling constant at each loop inte- 
gration according to Eq. O. The values of the form factor in the DL approximation of 



Ref. []I0| are also shown with their indicated scale of where was evaluated. Choosing 
an effective scale of roughly 9m^, one obtains results very close to the exact renormal- 
ization group improved values. 
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Figure 7: The effect of incorporating a running coupling constant at each loop inte- 
gration (circles) according to Eq. |l^ (suppressing the Sudakov term). The one- and 
two-loop running coupling solutions are within 8% of each other. Also shown are upper 
and lower limits according to the indicated values of in the DL approximation. 
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where a^Q is given in Eq. |T6| and asv is the exact one-loop resuh minus the one loop 
leading-logarithm pieces which are resummed in (t^q, i.e. 



o"sv 



<7SV,NLO — (^LO 



-6J^ + 



TT 



log 





(19) 



By adding the second (Sudakov) piece in the square brackets we remove (at least up to 
terms 0(/^/s) -C 1) the dependence on the gluon energy cutoff Ic- Note also that the 
complete expression for the two-jet cross section (with the remaining Ic dependence 
displayed) 

f^2j = (^RciQ + asY + (Th(/c, e, S) . (20) 

contains a mixture of exact 0{as) and resummed pieces. For the former, we use ml 
as the scale for a^-Q We believe that this is a more reasonable physical choice than the 
fixed value = 0.11 used for illustrative purposes in |11]]. The resummed contributions 
are based on the scale choice 1^ in the loops, as already discussed. 

Before computing and combining the various components of the two-jet cross sec- 
tion in Eq. ^ we must address the issue of the dependence on the unphysical infra-red 
parameter Ic- If we were to expand out the resummed RG-improved form factor o"|jq (/c) 
in powers of ^^(mg), and retain only the 0{as) term, we would find that the Ic depen- 
dence exactly canceled that of a}i{lc,e,6).^ However in the full resummed expression, 
there is nothing to cancel the explicit Ic dependence at higher-orders. The canceling 
terms would come from the as yet unknown higher-order contributions to (Th- Faced 
with this dilemma, we have several choices. We could, as in |]1T|, neglect the higher- 



order terms in the Sudakov form factor altogether, and include only the non-Sudakov 
form factor which is of course independent of Ic- Furthermore, we showed in [ITTI that 
with the choice e = (9(0.1), the combined contribution of virtual gluons and real gluons 
with Eg < e-\/s to was dominated by the non-Sudakov "6JF" part. This suggests 
that the most reasonable procedure for the resummed cross section is to take Ic ~ e^/i- 
We stress that this is an approximation, since it corresponds to making an assumption 
about the contribution of real multi-gluon emission with energies < 0{ey/s). 

As our 'best guess' RG-improved, resummed two-jet cross section, therefore, we 
have 

C^2j = (T^Si^y/^) + ^SV + C^H(ev^, 6, S) . (21) 

Figure ^ shows the y/s dependence of this cross section, normalized to the leading- 
order (Born) cross section. Evidently, the effect of the renormalization group improved 
results is significant when compared to the results in Ref. |jll| (where the scale of was 

^We chose the QCD scale parameter A such that as{ml) ~ 0.2 for mt — 4.5 GeV, at both leading 
and next-to-leading order. 
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This was shown explicitly in Ref. |11[|, see for example Fig. 3 therein. 
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a free parameter) . The Sudakov form factor reduces the all orders hard RG-form factor 
by over a factor of two, although it is still roughly a factor of three larger than if as{s) 
had been chosen in the hard DL-form factor. The hard Bremsstrahlung contribution is 
strongly dependent on the available phase space as the mass suppression is removed for 
these terms. The size of the Bremsstrahlung corrections is thus enhanced substantially 
through the large effective scale of the RG-Sudakov form factor. In total, we find that 
the RG effects lead to an enhancement of the background of roughly a factor of two 



compared with the scale chosen in Ref. |11] 



We note, however, that a similar effect for the Sudakov form factor is also present 
in the Signal H — > bb, so that the RG-effects in the non-Sudakov form factor-and in 
the misidentified cc background both Sudakov and non-Sudakov- will be decisive in a 
precise determination of the Signal/BG ratio. 



4 Conclusions 

In this paper we have shown how to incorporate the exact one and two loop running 
coupling into the massive Sudakov as well as the hard non-Sudakov double logarithmic 
form factors found in 77 (Jz = 0) — > qq. This process is very important at a future 
PLC as it contains the largest background to intermediate mass Higgs production. 
The RG-improvement is achieved without encountering the ominous Landau-pole sin- 
gularities in the hard form factor although the all orders leading and next to leading 
RG-solution is incorporated. The reason for this fortuitous fact is that the novel non- 
Sudakov form factors are "regulated" by the fermion mass, yielding an effective low 
energy cutoff. 

The size of the effect is significant and is found to stay within the theoretical upper 
and lower limits for both the form factor as well as the cross section contribution. The 
background is enhanced significantly compared to Ref. by the renormalization 



group effects as the effective scale is found to be roughly 9m^ for the hard form factor. 
The effective scale for the RG-Sudakov form factor is cutoff dependent. 

In this context it is worth pointing out a major difference to the 77 — > H signal 
contribution from a b-quark loop. There a function similar to the Abelian form factor 
of Eq. 1^ occurs, as was shown in Ref. |2^. There is, however, a big difference with 
respect to the renormalization of the strong coupling, as in the case of 77 — > H 
is only renormalized at the three-loop level. The effect of the running coupling in that 
process is thus much smaller, but perhaps still significant in the case of a large tanP 



enhancement scenario [21 



The effect of a running mass parameter is subleading as it occurs only as the 
arguments of logarithms. The mass term in the factorized Born amplitude is fixed 
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Figure 8: The total two-&-jet cross section (i.e. exact next-to-leading order contri- 
bution plus resummed Sudakov and non-Sudakov form factors) normalized to leading 
order, for the Sterman- Weinberg jet parameter choices (e = 0.1, 6 = 10° and 20°). Also 
shown is the separate non-Sudakov (dashed hne) and Sudakov (dash-dotted hue) form 
factor contributions, and their product (long-dashed line). 
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by the on-shell renormalization scheme employed in the exact one loop calculation of 
Refs. 

A full Monte Carlo study of both signal and background effects including inves- 
tigations of optimal experimental cuts and 6-tagging efficiencies is work in progress 
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5 Appendix 

In this appendix we give explicit results for the gluon mass dependent renormalization 
group improved massive virtual Sudakov form factor. We also give two loop results for 
both the virtual as well as the sum of the virtual plus soft real contribution and show 
numerically that both results exponentiate. 



5.1 The massive virtual RG-improved Sudakov form factor 

We begin with the result for J^Sv '- 



RG 



by 



ZTT JO Pi JO ai s s 
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1 + c log 

Cp { /"m d(i\ dai dj3i dai 



^ dPi dai dPi dai \ asim^) 

^ X%/3i^J 1 + clog^ 

f i log 4 flog ^ - lU i log ^1 (22) 



2tt [c \ as{s) J (? 

The A-dependent terms cancel out of any physical cross section when the real Bremsstrahlung 
contributions are added. In order to demonstrate that the result in Eq. ^ exponenti- 
ates, we calculate the explicit two loop renormalization group improved massive virtual 
Sudakov corrections, now containing a different "running scale" in each loop, and find: 
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s s J 

These integrals cannot be solved analytically anymore, however, numerically with the 
Monte Carlo generator Vegas it is straightforward. Fig. ^displays the numerical 
results of Eq. in comparison with half of the square of the RG- improved one loop 
results of Eq. The agreement is well within the statistical uncertainty and thus 
demonstrates that also the new running coupling contributions to the massive virtual 
Sudakov form factor exponentiates as expected. For completeness we also give the 
subleading terms of the pure one loop form factor which is of course also important 
for phenomenological applications. The complete result is thus given by: 

^RG a,{w?)Cp \\ s ( as{X^) ,\ , ^ . ^.(m^) 

1, s , 2tt''\ 

+ (24) 



5.2 Two loop results for the A-independent RG-Sudakov form 
factor 

We now give the two loop results for the physical combination J^g^ + 2JF|^ entering 
into cross sections. With the familiar Sudakov technique we find analogously to the 
purely virtual case: 
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s s 



18 



Virt. Exp. 



14 



13 



12 



11 



10 



O ; (lLf/2, 0.05 GeV^ 

A : 2L, 0.05 GeV^ 

□ ; (lLf/2, 0.1 GeV^ 

T ; 2L, 0.1 GeV^ 



m 



an 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



m 



100 120 140 160 180 200 



Figure 9: A comparison of the second term in the exponential obtained from the 
renormahzation group improved massive one-loop Sudakov form factor compared with 
the explicit two loop results. The quark mass is kept fixed at mj, = 4.5 GeV and the 
value for the gluon mass A is indicated in the figure. The result in Eq. ^ clearly 
exponentiates. 
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A numerical evaluation of Eq. ^5|-which contains a different "running scale" for each 
loop-is given in Fig. ^ and compared to half the square of the one loop result in Eq. ^. 
It can be seen that the renormalization group improved massive Sudakov form factor 
exponentiates as expected. 
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